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Abstract. - The linear threshold condition of the mirror mode including weak drift and finite Larmor radius 
effects is analysed and shown that it can be given to arbitrary precision in closed analytical form. Combined 
with observations of mirror modes in the magnetosheath it is shown that only magnetic 'dips' are produced by 
the mirror instability. A semi-experimental threshold condition is derived. The occasionally observed so-called 
'peaks' cannot be attributed to mirror modes but result probably from nonlinear evolution of the compressive 
(— I fast mode. 

O 

c/3 Introduction. - The magnetic mirror mode is one of the 
O^most interesting plasma modes in a collisionless anisotropic- 
O pressure high-/3 plasma. It is practically non-oscillating and 
c/5 occurs when the perpendicular (thermal) plasma (ion) pressure 
> sufficiently exceeds the parallel pressure. This mode 
'^has originally been inferred from simple fluid theory [1] and 
L_3n the following decades has been investigated more closely in 
ever higher precision (cf., e.g., [2-4]). It is of particular interest 
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because being of nearly zero frequency and being capable of 



assuming large amplitudes it may serve as energy source in the 
turbulent cascade that leads to low frequency magnetic plasma 
t — turbulence [5]. 

^ J Space plasma observations in various regions, planetary 

^SJ magnetosheaths, the inner terrestrial magnetosphere, near 
comets, and even in the solar wind provided evidence for the 
presence of mirror modes under the conditions prevalent in 

T\ these space plasmas. (While a comprehensive critical review 
^ of the observations and theory is in high demand, it is not the 

k>( place for it here. The reader is therefore referred to some less 
contemporary reviews, as e.g., [6].) Still, even though the mir- 
ror mode seems to be a very simple basic plasma eigenmode, 
controversy remains about its mechanism and physical realisa- 
tion. Several competing linear and nonlinear theories and some- 
times controversial opinions have been put forward. Numerical 
simulations have also been performed but are not in the posi- 
tion to clarify the case. It is, however, clear that for the mirror 
mode to occur as an eigenmode of anisotropic pressure plasmas 
it should evolve from a linear plasma instability starting from 
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the prescribed anisotropic initial conditions. These can be very 
different, either assimiing a homogeneous plasma, a stream- 
ing plasma, a weakly or strongly inhomogeneous plasma, a 
two-component plasma consisting of protons and electrons or a 
plasma with an admixture of heavy ions, and depending on the 
undisturbed ion distribution function and possibly even elec- 
tron temperature and electron pressure anisotropy. In addition, 
since the mode is a long wavelength mode being of macro- 
scopic size both along and perpendicular to the external mag- 
netic field, one might need to include the spatial boundaries of 
the region where the modes have been observed. This applies 
in particular to the Earth's magnetosheath which is bounded by 
the bow shock and the magnetopause. To account for all these 
initial conditions makes little sense. Hence one choses one or 
the other sufficiently simple initial condition. In order to clarify 
the basic state of the linear instabiUty it is most instructive to 
investigate the weakly inhomogeneous unbounded case in a bi- 
Maxwellian two-component plasma. The present Letter intends 
to clarify what can be learned from the available most elaborate 
linear mirror mode theory. We also have in view its application 
to regions in space where mirror modes have been observed but 
definitely restrict to the linear treatment only. Some conclusions 
can be drawn from its comparison with recent observations in 
the Earth's magnetosheath [7]. 

Instability threshold condition. - The most complete ki- 
netic derivation of the linear growth rate of the mirror mode 
including a weak inhomogeneity B(y) = Bq{\ — ay)z in an 
anisotropic high-j3 plasma under the condition that the plasma 
is bi-MaxweUian (at this stage of the theory an unimportant re- 
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striction) has been given by Pokhotelov et al. [4, 8-11 J. For 
the inhomogeneity to be weak only, the inhomogeneity pa- 
rameter is restricted as |a| <^ kj_, where k = (^x,^||) is the 
wave number of a plane wave. Since the inhomogeneity is 
strictly perpendicular to the magnetic field no such condition 
holds for . This contrasts other claims [12] of an apparent 
impossibihty to derive an instability threshold. In fact, k\i can 
become arbitrarily small. The smallness condition on the in- 
homogeneity also comes up for the decoupling of the mirror 
mode from the kinetic Alfven wave which on scales of the or- 
der of the ion inertial length A, = c/(Opi is another (non-zero 
frequency) eigenmode in high-j3 plasma. This coupling takes 
place at steep plasma gradients, i.e. it should play a role in the 
nonlinear evolution of the mirror mode whenever it grows to as- 
sume large amplitudes and evolves steep gradients. Since it has 
been shown that maximum growth of the mirror mode occurs 
at kj^pi = 0(1), with p,- = Vij^/cOd the ion gyroradius [8, 10] 
(v,x,v,|| are the respective perpendicular and parallel thermal 
velocities of the ions, and (Od is the ion cyclotron frequency), 
this condition also impUes that |ap, | <C 1. Linear theory then 
yields for the real frequency of the unstably excited mode 
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Here we introduced the pressure (or temperature) anisotropy 

A = (Tx/Tji) — 1 > which, for the mirror mode is always as- 
sumed to be positive. (Since for the mirror mode j3x is impor- 
tant, with J3|| being dynamically irrelevant, all expressions in 
this Letter are given in terms of A and jSx-) Under these as- 
sumptions the complete expression for the linear growth rate 
becomes 
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7o = (2/7r)z {|A;|||v,||/j3x(A+ 1)} is a positive factor. Clearly 
7 — > for k^ 0. This expression is independent of a. In the 
linear approximation the inhomogeneity enters only in the real 
frequency. 

InstabiUty arises for positive right hand side in Eq. (2). The 
first two terms, put to zero, yield the ordinary mirror thresh- 
old condition Aj3x > 1 for homogeneous plasma [2], neglect- 
ing the finite Leirmor radius effect in the third term. The fourth 
term gives a further correction obtained in the above reference. 
Both, the finite Larmor radius contribution and this term reduce 
the growth rate and increase the threshold. Clearly the clas- 
sical condition is an oversimplification of the real conditions 
even when neglecting the finite Larmor radius effect. Setting 
the right hand side in 7 to zero one obtains the threshold condi- 
tion for marginal instability of the mirror mode. 

Before drawing any conclusion about stability it is of interest 
to investigate the real frequency. This is clearly different from 
zero only in the case of inhomogeneity. Thus any modes have 
finite phase and group velocities in the plasma frame only when 
the inhomogeneity is taken into account or, otherwise, when the 
plasma is moving. However, Eq. (1) suggests that there are two 



possible modes whether Aj3x ^ 5 such that it seems to be nec- 
essary in the growth rate to consider these two cases separately. 
However, this is not the case as linearly growing waves in the 
plasma frame can be obtained only when the larger sign holds, 
as will become clear below. Growing drift-mirror modes prop- 
agate down the gradient in the direction of decreasing magnetic 
field, which intuitively is reasonable because of the decreasing 
magnetic pressure and stresses that support wave growth. 

Defining X = A + 1, it is easy to show that for the positive 
A > I the marginal stability threshold condition becomes 
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where we replaced J3|| = /3x / (A -|- 1 ) . The expression in the large 
bracket is positive definite. Hence, there is just one solution 
with X >0 which is given by 
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Here [• • • ] stands for the large bracket in Eq. (3). One immedi- 
ately realises that a slightly finite Larmor radius modified mir- 



ror threshold is obtained for small -C k]_ which reads 



Aj3x>l + ^feipf- 
with positive definite correction term 
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The larger J3x the larger is the effect of the term inside the 
bracket on the threshold, and for very large 1 the correc- 

tion term reduces to k^j2>k\. This case is unambiguously iden- 
tified as the Unear drift mirror mode. (The condition k^^ <C k\ is 
natural for linear mirror modes for which magnetic tension in- 
hibits oblateness. This holds even in the nonlinear case.) Both, 
the finite Larmor radius effect and the parallel wave number 
term act stabilising. Further improvements can be done by in- 
cluding finite electron temperature effects or anisotropic non- 
MaxwelUan distribution fiinctions, but do not change the con- 
clusion. 

Tuming to the case with A < j which describes uphill mov- 
ing waves, the threshold condition for instability becomes 
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This can be satisfied only when A < 0, a case that is not covered 
by the derivation of the mirror mode but applies to Alfven fire- 
hose waves. Hence, in this range mirror waves do not grow. Up- 
hill moving mirror modes are linearly stable. This does not pre- 
clude that they may not become parametrically excited when- 
ever the parametric growth rate exceeds the Unear damping rate 
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in Eqs. (2) and (7). This may happen when mirror waves evolve 
nonlinearly, generate sidebands and cascade into turbulence. 

Thus the sole case of linearly growing mirror modes is cov- 
ered by downhill propagating waves which, as is well known 
from linear analysis [4, 8, 9] take their free energy from res- 
onance with the bulk of low parallel energy ions of velocity 
V|| « 0. We note in passing that this bulk resonance makes 
the difference between mirror and electromagnetic cyclotron 
modes which otherwise grow under the same condition of pos- 
itive A > 0. The latter waves are in resonance with the high 
energy particle component of resonant energy S'res ^ mi{(Oci — 
a))^/2^J = - (o / (OciY / k^^Xf , where o) = a)(k) is the elec- 
tromagnetic ion cyclotron dispersion relation, A,- — c/(Opi the 
ion inertial length, and ~ niiV^/l the Alfven energy. Thus 
the two instabilities are about independent. They may or may 
not coexist. (Under low-j3 conditions the ranges of mirror and 
electromagnetic ion cyclotron instabilities have recently been 
investigated [15]). 

Application to magnetosheath observations. - In linear 
theory there is no difference between the mirror and drift- 
mirror instabilities in their growth rates. If streaming with 
velocity V would be included then the mirror mode would 
of course as well become oscillating with Doppler frequency 
0) = k • V [16] when transforming from the plasma to the labo- 
ratory frame. 

Observation of the mirror mode in its linear state is hardly 
possible. Wave spectra measured in the magnetosheath [16] 
have been interpreted as indicating the mirror dispersion rela- 
tion. However, the errors are very large and the linear approxi- 
mation is questionable as the observed waves have amplitudes 
b/B(, reaching from several % to several 10% indicating that 
one is dealing with nonlinear structures which develop steep 
spatial gradients in the magnetic field and plasma density and 
j3 , thereby violating the assumptions of linear theory and cou- 
pling the mirror mode to the kinetic Alfven wave. Observa- 
tions in the magnetosheath in addition indicated (e.g. [17]) that 
mirror structures seem to undergo an interesting spatial trans- 
formation between the bow shock, the outer boundary of the 
magnetosheath, and the magnetopause, the inner boundary of 
the magnetosheath. In the outer part of the magnetosheath the 
mode seems to nonlinearly develop into a sequence of 'mag- 
netic walls' (peaks), while closer to the magnetopause it occurs 
as 'magnetic holes' (dips). Admittedly, it is problematic to de- 
cide experimentally between these two forms simply because 
the definition of the background magnetic field in a highly 
structured (or spatially fluctuating) magnetic field is difficult 
to determine. 

Experimentally determined threshold. A recently per- 
formed more sophisticated and also more precise analysis [7] 
taking advantage of spatial Cluster four-point measurements 
in the magnetosheath confirms these observations. Soucek et 
al. [7] construct a quantity they call 'peakness' ^ with the help 
of which they can approximately distinguish between 'peaks' 
^ > and 'dips' 3^ <Q. This quantity correlates with jSy and 
jSxA (their Figure 2). If we eliminate from their data (which 
they provide as scatter plots) it is possible to construct the sta- 
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Fig. 1: The ranges of mirror 'dips' and 'peaks' in the plane of the pressure 
anisotropy A and /B^ for the Earth's magnetosheath based on the observations 
of [7] and purely linear theory. This representation makes use of the experi- 
mentally determined approximate boundary between 'dips' and 'peaks' (black 
line) and the theoretical mirror threshold condition. For comparison the crosses 
indicate the linear threshold condition for the miiTor mode. It is obvious that the 
min'or mode cannot produce any of the observed 'peaks' above the experimen- 
tal curve for anisotropics below A ~ 2.5. On the other hand, the experimental 
threshold allows for 'dips' both, above and below the linear mirror threshold. 

tistical relation jSy ~ 0.6 exp(0.3j3xA). Expressing jSy through 
the anisotropy this yields a relation between A and /3x which is 
of the form 

j3xA«0.5 + ln(^^^ (8) 

This expression is derived without distinguishing between 
'peaks' and 'dips'. It should be refined by considering these 
two data sets separately which might lead (in addition to a pre- 
cision of the numerical coefficients) to two separate relations 
of which only the one for 'dips' will apply to the mirror mode. 
This will become clear below. 

In spite of these reservations, the above relation is interesting 
as it resembles a 'nonlinear threshold condition' for the mirror 
mode in close similarity to the simple linear threshold condi- 
tion. The logarithmic nonlinearity in this expression is a con- 
sequence of the notoric observation of mirror modes in their 
nonlinear state. However, the nonlinearity is logarithmic and is 
thus weak. The main difference is in the statistically uncertain 
numerical term 0.5 on the right. (Unfortunately, in [7] no lower 
limit on I ^1 has been set such that the correlations are infected 
by the accumulation of irrelevant cases around ^ « 0.) 

'Peaks' or 'dips'? Only 'dips' survive. Analysing their 
distribution of 'dips' and 'peaks', the authors of [7] found 
an approximate relation 7x/7]| « 2.2/^^-^ between the tem- 
perature anisotropy and jSy for the boundary between the re- 
gions of 'dips' and 'peaks' [7]. For values below this approxi- 
mate boundary curve they predominantly observe 'dips', while 
above it they predominantly identify 'peaks'. Rewritten in our 
terms for the 'dips' this relation becomes 

A< (3.5/j3]'''4)-l (9) 

The numbers in this expression (in particular the power on j3x) 
are quite uncertain, however, due to the large scatter in the data 
[7]. (Note that the power of j3x is very close to |.) Still, this 
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can be compared with the marginal mirror instabihty threshold, 
yielding for 'dips' 

1 < AJ3^ < 3.5 /3 0-36 - /3i ~ 3.5 jS] - J3^ (10) 

This condition can be satisfied as long as 1 < J3_l < 6.3, which 
yields for the anisotropy that 0.15 < A < 2.5 for the occur- 
rence of 'dips'. Thus both, anisotropy A and ji±, are restricted 
to small values. 

On the other hand, accounting for 'peaks', the above inequal- 
ity Eq. (9) is inverted, and all anisotropics A > 2.5 are permit- 
ted for causing 'peaks' . These regions are shown graphically in 
Figure 1 where we also have drawn the linear mirror threshold 
boundary (crosses). 

Given that the experimentally determined condition [7] 
holds, we apparently find that it is in accord with the approx- 
imate threshold relation of the mirror mode criterion for mod- 
erately large j3j_ to cause 'dips', while the peaks observed by 
Soucek et al. [7] to accumulate in the large white region above 
the boundary in Figure 1 obviously are no mirror mode signa- 
tures. Only peaks at large anisotropies could possibly be caused 
by mirror modes. However, from linear analysis this can hardly 
be definitively concluded to be the case. It is also interesting to 
see that the allowed region of 'dip' occurrence centres around 
the marginal threshold condition of the linear mirror mode. It 
is suggestive that 'dips' are found both above threshold (in the 
linearly unstable domain) and below threshold (in the linearly 
stable domain), however with both regions being quite narrow 
in parameter space. 

Discussion. - In order to reach the observed large 'dip' 
ampUtudes the mirror instabihty needs to enter the nonhnear 
state, which might have proceeded along the lines investigated 
for evolution near marginal stabihty [18, 19]. The main con- 
clusion we can draw from linear theory is that it is improbable 
that the mirror instability near threshold would be responsible 
for generating 'peaks' . This conjecture could experimentally be 
proven by eliminating all statistically uncertain small amphtude 
I ^ I values from the data of [7] and investigating the data sets 
for 'peaks' and 'dips' separately. 

On the other hand, nonhnear mirror mode theory [19] and 
simulations [20, 21] give a clue on the evolution of the mir- 
ror instability. It has been proposed recently in [18] that near 
marginal stability mirror mode collapse could take place. If this 
happened it would completely determine the further evolution 
of the instability, in particular when it starts from a large am- 
plitude disturbance. Though this theory is very attractive and 
beautiful, a more sophisticated and elaborated nonlinear inves- 
tigation [19] has shown that the terms in the nonlinear evolution 
of the mirror instability that lead to collapse do in fact cancel 
such that collapse does not take place. Instead, the mirror mode 
evolves solely into 'dips', in perfect agreement with the results 
of this investigation which are based exclusively on hnear the- 
ory combined with observation [7]. Saturation of the mirror 
mode in its nonlinear state is then provided by trapping of the 
resonant low-parallel energy ion component [19] that causes 
the magnetic field to flatten out in the centre of the mirror 



'dips' , similar to what had been found in the very high-time res- 
olution magnetic observations of mirror modes [22] provided 
by the Equator-S spacecraft. 

What concerns the 'peaks', it seems that from a theoreti- 
cal point of view the mirror instabihty is not responsible for 
their evolution, neither when it evolves from thermal fluctua- 
tion level, nor when it starts from large amphtude disturbances. 
In the latter case the mirror mode becomes the drift mirror 
mode and preferentially propagates into the direction of de- 
creasing magnetic field thereby further depressing the field. 
Uphill moving modes could possibly be generated by paramet- 
ric instability requiring the existence of large amplitude 'dips'. 
In this case they would be short wavelength structures such as 
have been found in the Equator-S recordings inside mirror bub- 
bles. Compressions of the field should be expected to occur 
in contrast when a compressive wave mode is excited in the 
plasma. This is naturally the case when the fast mode is unsta- 
ble, a case that was suggested also in [19]. The observation that 
compressive 'peaks' in the magnetosheath are found close to 
the bow shock hints on a connection to the shock as, in particu- 
lar, at short distances from the shock the thermal anisotropy of 
the magnetosheath plasma has not yet built up sufficiently and 
the plasma is probably mirror stable as was suggested in [17]. 
From quasi-parallel shock theory one may indeed expect that 
behind the quasi-parallel bow shock the shock radiates com- 
pressive waves which in the magnetosheath appear as large am- 
phtude compressions or sohtary structures on the fast mode. 

* * * 
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